Role of electronic nematicity in the interplay between s- and d-wave broken-symmetry 
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To understand the role of electronic nematic order in the interplay between s- and d-wave particle- 
particle or particle-hole condensate states, relations between various s- and d-wave order parameters 
are studied. We find that the nematic operator transforms two independent six-dimensional vectors. 
The d-wave superconducting, d-density wave, and antiferromagnetic orders are organized into one 
vector, and the s-wave superconducting, charge density wave, and spin-triplet d-density wave orders 
into the other vector. Each vector acts as a superspin and transforms under the action of SO(6) 
where charge, spin, r\- and 7r-pairing, spin-triplet nematic operators satisfy the SO(6) Lie algebra. 
Electronic nematic order is not a part of the SO(6) group. It commutes with all 15 generators. Our 
findings imply that nematic order does not affect the competition among the order parameters within 
the same superspin, while it strongly interferes the interplay between two order parameters that 
belong to different superspins. For example, nematicity allows a linear coupling between d- and s- 
wave superconducting order parameters which modifies the superconducting transition temperature. 
A generalized Ginzburg-Landau theory and further physical implications are discussed. 

PACS numbers: 



INTRODUCTION 

A minimal model for strongly correlated materials in- 
cluding the most complex systems such as the high tem- 
perature cuprates is the Hubbard or t-J model. However, 
even within the simplified Hubbard model away from the 
half-filling, there has been no consensus about the ground 
states, and various theoretical proposals have been made 
for the phase diagram of the high temperature cuprates. 

One of such states, the staggered flux phase (d-density 
wave), was discussed in exitonic condensation [1], t-J 
model, and Hubbard model [2-5], and further suggested 
as a pseudogap phase of the cuprates [6]. The spin-triplet 
version of the staggered flux phase (spin-triplet d-density 
wave) was also discussed in the context of high tempera- 
ture cuprates. [7-9]. 

Another proposed particle-hole condensate state of the 
angular momentum 1=2 channel with broken rotational 
symmetry is the electronic nematic state. [10-18]. The 
spontaneous formation of nematicity has been invoked 
to explain the anisotropic transport observed in a two- 
dimensional electron gas in a high magnetic field [19, 20] 
and in Ru-based oxides [21]. Its relevance to high tem- 
perature cuprates was evidenced by a neutron scatter- 
ing measurement where anisotropic scattering patterns 
have been observed in YBa2Cu306.5.[22] From a weak- 
coupling point of view, it is sometimes also called Pomer- 
anchuk instability [23]. 

Given various proposed order parameters, the inter- 
play between s- and d-wave order parameters has been 
of intensive theoretical study. Examples for s- and d- 
wave order parameters are s- and d-wave superconduc- 
tors, charge density wave, spin density wave, d-density 
wave, spin-triplet d-density wave, and nematic states. 



Among them, it was reported that nematicity plays an 
important role in the interplay between s- and d-wave 
superconductors [24], and the spin density wave and spin- 
triplet d-density wave states [9]. However, a full set of 
the relations between them is still lacking. 

In this paper, we offer a complete theory on how s- and 
d-wave orders transform via the nematic order, and rela- 
tions between them. We found that the order parameters 
listed above can be organized in two independent six- 
dimensional vectors. One vector is composed of d-wave 
superconducting, d-density wave, and spin density wave 
order parameters, while the other vector contains s-wave 
superconducting, charge density wave, and spin-triplet d- 
density wave order parameters. Each vector transforms 
under the action of SO (6). Charge, spin, //-pairing [25], 
7r-pairing[26], and spin-triplet nematic [27] operators to- 
gether satisfy the SO (6) Lie algebra. 

The nematic order parameter commutes with all gen- 
erators, and hence is not a part of the SO (6) group. How- 
ever, it transforms the two independent vectors connect- 
ing s- and d-wave order parameters. Our findings imply 
that nematic order does not interfere the competition be- 
tween the order parameters within the same vector, but 
strongly affects the interplay between two order parame- 
ters belonging to different vectors. For example, it allows 
a linear coupling between s- and d-wave order parameters 
in the Ginzburg-Landau (GL) free energy, which modi- 
fies the physical properties of both phases. We found 
that the conditions for non-zero linear coupling differ for 
particle-particle and particle-hole condensate states. 

Below we will review an SO (6) group theory and 
present the relations between the nematic order, the gen- 
erators, and superspins. We will show the role of nematic 
order for particle-particle condensate states, and elabo- 
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rate a similar process for particle-hole cases. We will also 
discuss the implications of such relations in the context 
of GL free energy, and superconducting transition tem- 
perature related to high temperature cuprates. 



NEMATIC ORDER PARAMETER AND SO(6) 
GROUP 



It was first pointed out by Yang[25] that the 
77-pairing state is an eigenstate of the Hubbard 
model. The 77-pairing operator is defined as i) + = 

-«Ekaa' c k,<7°'L' c -k+9, ( 7' and V~ = where Q = 

(71", 7r) and a y is a Pauli matrix. The 77 operator carries 
charge 2e and spin 0, and commutes with the Hubbard 
Hamiltonian at half filling fi = U /2, where \i is the chem- 
ical potential and U is the on-site Hubbard interaction. 
It is also an eigenstate of the momentum operator with 
the eigenvalue Q. 

Later it was found that the 77-pairing operators com- 
bined with the charge operator satisfy an SU(2) alge- 
bra (named pseudospin)[28], and further recognized that 
the Hubbard model has two sets of commuting SU(2) 
symmetries. One set is characterized by the pseudospin 
of 77-pairing and charge operators, and the other is con- 
ventional spin operator. [29] The three-dimensional vector 
transforming under the action of pseudospin SU(2) forms 
a superspin, where its three components are s-wave su- 
perconductor, and charge density wave with the ordering 
wave vector Q. It was also reported that the pseudospin 
SU(2) rotates another superspin composed of d-wave su- 
perconducting and d-density wave order parameters. [30] 

On the other hand, the vector transforming under the 
spin SU(2) is the spin density wave with the ordering 
wave vector Q. Under a particle- hole transformation for 



operators: 



one spin species, q, — > (— lfc^, the role of the two sets 



of SU(2) generators is interchanged. The same particle- 
hole transformation maps the positive Hubbard model 
to the negative Hubbard model, and it also maps the 77- 
pairing to the Nagaoka ferromagnetic state. [31] A further 
generalization of the concept of the exact SO (4) symme- 
try of the Hubbard model to a unified theory of anti- 
ferromagnetism and d-wave superconductivity based on 
SO (5) symmetry was later proposed to understand the 
physics of the high temperature cuprates. [26, 32] 

Here we present a full list of relations between s- and d- 
wave order parameters including those mentioned above. 
The ground state order parameters transform as a six- 
dimensional vector under the action of SO (6). They can 
be organized into a vector n a (a — 1...6) which should 
satisfy 

[Lab, n c ] = -i (S bc n a - 5 ac h b ) . (1) 
where L a b are generators of SO(6) based on the following 
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where k runs over the reduced Brillouin zone, d(k) = 
cos k x — cos fcy, a takes the values x, y, z, and a a are the 
Pauli matrices. 

The generators of SO (6) can be represented by an an- 
tisymmetric 6x6 matrix, L ab = —L oa . 
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where 5ft C> = \(6~ +6+) and SO = 
satisfies the correct SO(6) Lie algebra, [33] 



6+). It 



i ( SadLbc + SbcLad — SbdLac — S ac L 



bd 

Here L12 is the charge operator, and £34,^35, and L45 are 
the three components of the spin operator. Liq and L26 
represent real and imaginary part of the 77-pairing. [25] 
Li3, L14, and L15 (L23, £24 and L25) denote x, y and z 
component of the real (imaginary) part of the 7r-pairing 
which carries charge 2e and spin 1 and represents a bro- 
ken translational symmetry. [26] £36, L46, and L 5e corre- 
spond to the spin-triplet nematic order parameter carry- 
ing spin 1 and representing a broken x-y symmetry on a 
square lattice. [27] 

There exist two independent vectors. Each vector 
acts as a superspin, and transforms under the SO (6). 
This observation was first reported in Ref. [34], and 
a similar SO (6) symmetry was found in Fe-pnictide 
superconductors [35]. One superspin (superspin-1) con- 
sists of spin-density wave (A sa < w ), d-density wave (A^^), 
and d-wave superconducting (Aj; sc ) order parameters: 
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where hi = 0? A dsc , h 2 = 3 A dsc , h 3 = A* dw , h 4 = A v sdw , 
h 5 = A z sdwl and h 6 = A ddw . 

The other superspin (superspin-2) rotated by the same 
15 generators L a {, is composed of spin-triplet d-density 
wave (Atsf), charge density wave (A cdw ), and s-wave su- 
perconducting (A ssc ) order parameters: 



A " s / = 2 E d ( k ) RX^+Q,*' 

kaa' 



C k+Q,cr' J <T<T' Ck .<'' 



,+ 
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. ^ C ik,t C -k,i + C ik+Q,t C -k+Q4 

k 



(A+ ) f , 



Acdu, — ( c k 



: J<r Ck+Q,a + c ik+Q,<T C k,<T 



(5) 



where the superspin is arranged as hi = 5RA SSC , hi = 
SA SSC , h 3 = Af sf , h 4 = A v uf , h 5 = A z sf , and n 6 = 

A cdw . 

What is the role of the nematic order parameter within 
the SO (6) group? The nematic order parameter is given 
by 



N = E d ( k ) ( C k, CT C k,<x - 4+Q l( r C k+Q^ 



(6) 



When (N) = A/q 7^ 0, the phase is characterized by a 
broken x-y symmetry of the square lattice, and it is triv- 
ial to generalize to a broken point group symmetry in 
other lattices. Note that the nematic order parameter 
commutes with all 15 generators: 



[N,L ab } = 0. 



(7) 



This means that nematicity is not an SO (6) symmetry 
breaking field, and does not interfere with the compe- 
tition between the order parameters within the super- 
spin. For example, the phase diagram between antiferro- 
magnetism and d-wave superconductivity (both belong 
to superspin-1) studied in the t-J model based on SO(5) 
symmetry [26] (a subset of the SO (6) in this study) is not 
modified by the presence of nematicity. 

However, the nematic operator does not commute with 
the following conventional quantum rotor model. 
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q r -^Y. L\ab + E 



r a hih^, 



(8) 



i,a<b 



where the first term is the kinetic term and \ is the 
moment of intertia, and the second term is the poten- 
tial term. The Hamiltonian has SO(6) symmetry when 
r a is idential to all a. Note that the nematic operator 
commutes with the first term, but not the second term. 
On the other hand, the competition between the order 
parameters in the same superspin h a is determined by 
difference in r„ . 



What are relations between nematic order and other 
order parameters? The nematic operator transforms the 
components of the two independent superspins as follows: 



^tsf ' ^sdw 



= 2N, 
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= 2N, 
= iN, 



A~ A+ 



\N. (9) 



The nematic operator transforms s- to d-wave order pa- 
rameters which belong to different superspins. 

The above results are summarized in the table below. 



SO (6) generators 


Q, S, i)-, 7r-pairing, spin nematic operators 


nematic operator 


commutes with generators & 




transforms superspin-1 and -2 


superspin-1 


dSC, dDW, SDW 


superspin-2 


sSC, CDW, spin-triplet dDW 



TABLE I: A summary of the SO (6) group and the relations 
to nematic order. 



In the following section, we discuss the physical im- 
plications of the commutation relations using a GL free 
energy theory assuming that nematic order is present. [36] 



GINZBURG LANDAU THEORY 

The commutation relations in Eq. 9, [A, B] = N, in- 
dicate that if (N) = Ao is finite, a linear coupling be- 
tween A and B phases, such as 7 $ W with $ = (A) and 
\& = (B) , may be present in the GL free energy. The GL 
free energy is then given by 



T 



-$ 2 + 7 * $ + 



v$ 4 



(10) 



Assuming that a > 0, b > 0, and ab > j 2 (none of the 
phases represented by $ and is ordered), the solutions 
of the two coupled equations for $ and \& leads to the 
following dispersion of modes: [37] 



X ^(k) = £l(k) + £2(k) ± iVMk)- £2 (k)) 2+ 47 2 



where 



d(k) = a + p{(l + Af )k 2 + (1 
e 2 (k) = b + p{(l+Af )k 2 x + (l- 



-M))ky}, 

-Mo)k 2 }. 



(11) 



(12) 



Here we have used the effective Lagrangian L e ff = 
f (d t $) 2 - § {(1 +A/-0XW + (1 - A/o)(c\$) 2 } ~ f * 2 + 
§ (d t ^) 2 - §{(1 +AT Q )(d a y) 2 + (1-Af )(d y ^) 2 } - |* 2 - 
7^$, where p is the stiffness. Note that the excitations 
are anisotropic due to nematicity [38], and k x and k y are 
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deviations from an ordering wave-vector which is either 
or Q depending on the nature of ^ (or $). 

One of the solutions becomes when 7 = >/ab, lead- 
ing to an ordered state. The condensed state is a linear 
combination of VP and $, and the dominant contribution 
depends on a and b. Also, if one of them, say *5f, is finite 
(when a < 0), the other, $, is always induced as long as 
7 is finite. 

Is 7 always finite if nematic order exists? For exam- 
ple, consider a system in the nematic state with SO(6) 
symmetry at high temperatures. At low energy, the sys- 
tem spontaneously breaks the SO(6) symmetry, and one 
of the phases represented by "J is stabilized. If ^ repre- 
sents the d-wave superconducting state, $ is the s-wave 
component. Similarly if \1/ is the spin density wave, $ 
should be the spin triplet d-density wave. [9] Does ne- 
maticity always lead to an induced order parameter of 
$ without any extra condition? Below we show that it 
requires another condition for a non-zero linear coupling 
(in addition to the nematic order), and that the condition 
for a finite 7 differs for particle-particle and particle-hole 
condensates. 



DIFFERENCE BETWEEN 
PARTICLE-PARTICLE AND PARTICLE-HOLE 
PAIRS 

Let us compute 7 for particle-particle condensate 
states. To check the condition for a non-zero linear cou- 
pling coefficient 7 between d-wave and s-wave supercon- 
ducting cases (\f r = Re(A c [ sc ) and $ = Re(A ssc )), wc 
introduce ip^ = (c ko ., c_k-<x)- Then the order parame- 
ter is written as A ssc = J^k W^i^k- Inside the nematic 
state, the quasiparticle Green's function is written as 



G : (k, iu) n ) = -ioj n + 



(13) 



where 



£k = — 2£(cos k x + cos k y ) + 2id(k)Ao — At' cos k x cos k y , 

and (i is the chemical potential, t and t' represent the 
nearest neighbor and second nearest neighbor hoppings, 
respectively. Assuming that d- and s-wave superconduct- 
ing fluctuations couple to fermions with interactions of g\ 
and 52, the 7 coefficient becomes 



Ids 



5l52 T^d(k)^Tr(G(k, 



»g(£fc) - n F (-£ k ) 
2^ 



(14) 



where £ k = e k — [i. 7 is always finite as long as /x and/or t' 
is finite. In other words, when the particle-hole symmetry 
is broken and nematic order is present, the linear coupling 



term induces d- or s-wave superconducting order as we 
discussed in Eq. 10. 

However, the above result is not true for particle-hole 
pairs. 7 then is zero independent of particle-hole sym- 
metry. To examine the condition for particle-hole cases, 
let us introduce i/> kcr = (c k(T , c k+ Q CT ). In this basis, the 
Green's function becomes 

G _1 (k,iw„) = - ioj n I + e k T 3 - /i k /, (15) 

where e k = — 2t(cosk x + cosk y ) + 2<d(k)Ao = — ?k+Q 
and /ik = At 1 cos k x cos k y + fx = /^k+Q • The 7 coefficient 
for example between the charge density wave and the 
d-density wave order is then obtained as 

'Jcdw — ddw C< 

d(k) ^2 Tr (G(k, zu; n )TiG(k, iio n )r 2 ) 

k iuj n 

= 0. (16) 

This is similar to the coupling between ReAd sc and 
ImA ssc . This linear coupling is not allowed in the free 
energy due to symmetry. 

Let us consider the coupling between different direc- 
tions of spin density wave and spin-triplet d-density 
wave. For example, the coupling between antiferromag- 
netic fluctuations along the x-direction and spin-triplet 
d-density wave fluctuations along the y-direction are 
given by 8A x sdw oc ip\ji^ and SA^ sf oc Vk r iV'k where 

V'k = ( C ikt' c k+Qi)- Then the coefficient 7 is obtained 



Jtsf-sdw oc d(k)Tr (riG(k»o; n )TiG(kia; n )) 



E (n F (ik - /Uk) - n F (-e k - /Zk)) 
ze k 



0. 



(17) 



Note that the coupling is also 0, because both d(k) and 
e k change sign under k — > k + Q, while ^ k does not. The 
physical reason is that the spin density wave state breaks 
time reversal symmetry, while the triplet staggered flux 
does not. This fact is reflected in the commutation rela- 
tions, where Eq. 9 has the imaginary factor i. Therefore, 
a linear coupling is not allowed between s- and d-wave 
particle-hole condensate states. 

However, in the presence of a magnetic field h, the re- 
sult alters. Note that the particle-particle and particle- 
hole order parameters are related by the particle-hole 
transformation, which also maps the chemical potential 
to the magnetic field to be discussed in detail below. In 
the presence of an external magnetic field, 7 changes to 



ltsf-sdw(h ^ 0) 
d(k) 



(18) 



E 



2(e k + h) 



(n F (e k + h-/j,)- n F (-e k -h-fx)) 
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7 between the spin-triplet d-density wave and spin den- 
sity wave is finite when Ao and h are finite. The leading 
contribution of h and Ao to 7(Ao,/i) can be written as 
7(A/"o, h) = 70A0/1, where 70 depends on the interactions 
between fermions and the fluctuations of the order pa- 
rameters. [9] 

To understand the difference between particle-particle 
and particle-hole condensates, let us consider the 
particle-hole transformation. The particle-hole transfor- 
mation mapping the positive to the negative Hubbard 
model discussed above maps each component of the su- 
perspins as follows: 

-» ( A L^»). (19) 

where a = x,y, z. 

In addition to the known result that the antiferromag- 
netic order transforms to the s-wave superconducting and 
charge density wave orders, we found that the spin-triplet 
d-density wave phase transforms to the d-wave super- 
conducting and d-density wave orders, while nematic or- 
der is invariant. Since the chemical potential maps to 
the magnetic field under the particle-hole transformation, 
the conditions for a finite 7 between particle-particle and 
particle-hole condensates are also related by the particle- 
hole transformation - 7oA r /Lt(5RAd sc )(3ff A ssc ) maps to 
7oA r /i(A^ s ^)(Ag du; ) under the particle-hole transforma- 
tion. Therefore, the linear coupling between d- and 
s-wave superconducting order parameters requires a fi- 
nite chemical potential, while the coupling between spin- 
triplet d-density wave and spin density wave requires a 
magnetic field. Note that both breaks SO (6) symmetry, 
as the chemical potential and magnetic field appear as 
/1L12 and /1-L34 in Hamiltonian, respectively. 

EFFECT ON SUPERCONDUCTING 
TRANSITION TEMPERATURE 

Let us reexamine the GL free energy, Eq. 10, to see 
if the superconducting transition temperature is modi- 
fied by the coupling between the d- and s-wave super- 
conducting order parameters. We consider * = ($lA dsc ) 
and $ = (3fA ssc ), and 7 is finite and proportional to 
the nematic strength Ao , and particle- hole symmetry is 
assumed to be broken. 

Since the chemical potential couples to the charge oper- 
ator L12, it favors the d-wave superconducting state over 
the antiferromangetic and d-density wave states. The 
competition between antiferromagnetism, d-wave super- 
conductor, and d-density wave is determined by SO(6) 
symmetry breaking terms, where nematicity does not af- 
fect the interplay between them. 

Assuming that the superconducting state is stabilized 
in a finite window of phase space, and the transition tem- 



perature is set by T® ( a < below T° and assume b > 0), 
we are interested in the effect of nematicity on the super- 
conducting transition temperature. It is straightforward 
to check that the effective mass term a e ff is modified 
by a - ^ after integrating out the $ field. Note that 
a oc (T — T®) and a e j f oc (T — T c ) where T c is the transi- 
tion temperature modified by the coupling 7. Since the 
effective mass gets smaller due to the coupling to the s- 
wave component, the transition temperature T c is higher 
than T°. However, one should note that the current de- 
scription is based on a classical theory, and quantum fluc- 
tuations beyond the present study should be taken into 
account to see if the result may qualitatively change. 

DISCUSSION AND SUMMARY 

We have studied the role of the nematic order pa- 
rameter in the interplay between s- and d-wave particle- 
particle or particle-hole condensate states. These con- 
densate states include d- and s-wave superconductors, 
d-density wave, spin-triplet d-density wave, spin density 
wave, and charge density wave phases. We found that the 
nematic operator transforms d- to s-wave superconduc- 
tors, spin-triplet d-density wave to (s-wave) spin-density 
wave, and d-density wave to (s-wave) charge-density wave 
operators. This can be summarized as a transforma- 
tion between two different six-dimensional vectors. One 
vector is composed of d-wave superconductor, d-density 
wave, and spin-density wave order parameters, while the 
other vector consists of s-wave superconductor, charge- 
density wave, and spin-triplet d-density wave order pa- 
rameters. Each vector acts as a superspin and transforms 
under the action of SO(6). There exist 15 generators, 
which correspond to charge, spin, spin-triplet nematic, 
rj- and 7r-pairing operators, which form the SO(6) group. 

The transformation between the two superspins via ne- 
maticity implies that a linear coupling between two or- 
der parameters that belong to two different vectors can 
be present in the GL free energy. Such a linear coupling 
allows induced ordering when one of them is condensed. 
However, we found that there is an additional condition 
for a non-zero linear coupling, which differs for particle- 
particle and particle-hole condensates. For example, 
when d-wave superconductor (particle-particle conden- 
sate) and nematic order coexist, s-wave superconducting 
order is induced, only when the particle-hole symmetry 
is broken. On the other hand, when spin-density wave 
(particle-hole condensate) and nematic order coexist, a 
similar transformation allows an induced spin-triplet d- 
density wave, only when time-reversal symmetry is bro- 
ken. These results are consistent with symmetry consid- 
erations. Since the spin-triplet d-density wave does not 
break time reversal symmetry, while spin-density wave 
does, a linear coupling between the two order parame- 
ters is allowed when time reversal symmetry is broken by 
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an external magnetic field. 

It is also interesting to notice that the nematic operator 
commutes with the generators. When the Hamiltonian 
contains a term —g^2ij NiNj which favors nematic or- 
dering, it does not act as an SO(6) symmetry breaking 
field. It means that the nematic order can exist without 
interfering the competition among the six different order 
parameters within a superspin. It is merely a spectator. 
However, it affects the interplay between order parame- 
ters which belong to two different superspins. Nematicity 
allows a linear coupling between the two order parame- 
ters, and affects the physical properties of both phases. 
As an example, we showed that the d-wave superconduct- 
ing transition temperature is modified by the coupling to 
the s-wave superconducting order parameters which hap- 
pens when nematicity is present and particle-hole sym- 
metry is broken. 

The nematic order parameter has been widely dis- 
cussed in the context of strongly correlated materials. 
In particular, the phase diagram of the high temperature 
cuprates is complex and its complete understanding re- 
quires further experimental and theoretical investigation. 
Our results indicate that the proposed nematic phase af- 
fects phenomena in the superconducting phase such as an 
anisotropy in the spin susceptibility and an increase in 
superconducting transition temperature. It also affects 
antiferromagnetism via the coupling to the spin-triplet 
d-density wave when a magnetic field is applied. We 
do not attempt to find a microscopic Hamiltonain with 
SO(6) symmetry which is beyond the scope of the current 
study. However, we emphasize that Eq. 7 and 9 are ex- 
act independent of symmetry of Hamiltonian, and SO(6) 
symmetry is useful to identify the compact relations be- 
tween the nematic and other order parameters suggested 
in the context of high temperature cuprates. The GL 
free energy analysis hints the importance of nematicity 
for the phase diagram of antiferromagnetism and d-wave 
superconducting phase. 
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